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Abstract. We have detected relativistic and/or kinematic post-Keplerian 
parameters of no fewer than five pulsar-white dwarf binary systems in 
long-term timing observations at Arecibo. We discuss the resulting con- 
straints on pulsar and companion masses. The general trend is that the 
pulsar mass measurements are compatible with the canonical value of 
1.35 Mq, but in most cases somewhat higher masses are preferred. 



1. Introduction 

Evolutionary considerations suggest that pulsars in white dwarf-neutron star 
binaries might have masses higher than the canonical value of 1.35 Mq, perhaps 
due to extended, stable mass transfer. Here we report on observations of five 
pulsars in such binary systems. We have measured one or more post-Keplerian 
effects in each system, leading to constraints on the pulsar masses. (For a review 
of previous pulsar mass measurements, see Thorsett & Chakrabarty 1999.) 

We observed the pulsars over several years with the post-upgrade Arecibo 
telescope. We used the Princeton Mark IV data acquisition system to perform 
baseband sampling, coherent dedispersion, pulse-synchronous folding, and tim- 
ing (Stairs et al. 2000). Depending on the pulsar, we also used data collected 
at Green Bank (140 Foot telescope), Jodrell Bank, and/or Effelsberg, as well 
as data from the Arecibo Berkeley Pulsar Processor and the Princeton Mark III 
system (pre- upgrade Arecibo). Details will be given elsewhere. 



2. Constraints on Pulsar and Companion Mass 
2.1. Keplerian Mass Function 

For each system, five Keplerian parameters are measured: orbital period. Pi,] ec- 
centricity, e; angle of periastron, oj; time of periastron passage, Tq; and projected 
semi-major axis, x = (ai sini)/c, where oi is the semi-major axis of the pulsar 
orbit, i is the orbital inclination angle, and c is the speed of light. These com- 
bine to form the mass function, fi = (m2 sinz)'^(mi + m2)~^ = {2tt/ Pb)'^x'^ /Tq, 
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where Tq = 4.925 fis is the solar mass in time units, and where mi and 772-2 
are the mass of the pulsar and its companion, respectively, in solar mass units. 

While i is unknown, the constraint sinz < 1 leads to nii < rrul'^ — m^. 

2.2. Relativistic Phenomena 

Three relativistic phenomena constrain the system masses. 

Precession of an orbit can be measured if eccentricity is sufficiently large. 
(For circular orbits, precession is completely covariant with orbital period.) The 

rate of apsidal advance is cj = 3(P5/27r)^^/'^(l — e^)~"'^TQ^(mi+m2)^'''^, so mea- 
surement of Co yields total mass, mi + m2. 

Shapiro delay, the excess propagation time through the gravitational poten- 
tial well of the secondary is, for small e, Ats = — 2m2TQ ln[l — sinzsin((/) — 0o)], 
where (j) is the orbital phase measured from the ascending node, (f)Q. For edge- 
on orbits (large i), detection of Shapiro delay leads to measurement of m2 and 
sini. For face-on orbits (small i), Shapiro delay is covariant with projected semi- 
major axis and cannot be detected; this non-detection can be used to exclude 
high-inclination, low-companion-mass orbits. 

Orbital decay due to gravitational radiation causes Ph to change at a rate 

n = -(192V5)(Pfe/27r)-5/3(l+I|e2+§e4)(l-e2)-7/2T(^/Vim2 (mi+ma)-^/^. 
This can be measured for sufficiently small orbits (low Pb). 

2.3. Kinematic x 

Motion of a binary system relative to the Sun induces a change in the inclina- 
tion angle of the orbit, and hence a change in the projected semi-major axis, 
x/x = — /i cot i sm.6, where is the (measured) proper motion and 6 is the 
(a priori unknown) difference between the position angle of proper motion and 
the position angle of the ascending node (Kopeikin 1996; see also Nice, Splaver, 
& Stairs 2001). Measurement of x alone does not yield a value of i, but the 
constraint | sin0| < 1 leads to i < tan~^(^x/i;), giving a firm upper limit on i. 

2.4. Orbital Period Core Mass relation. 

Studies of low- and intermediate-mass binary evolution predict a nearly unique 
relation between and m2 in pulsar-helium white dwarf binaries (Podsiad- 
lowski, Rappaport, & Pfahl 2002; Tauris & Savonije 1999; Rappaport et al. 
1995). The Pb — m2 relation may be used to constrain the system masses; alter- 
natively, measured system masses may be used to test the evolution theory. 

3. Timing Analysis and Results 

Table 1 summarizes the post-Keplerian phenomena detected in each of five pulsar 
binaries. For each pulsar, we used tempoEI to calculate the goodness-of-fit of 
timing solutions over a uniform grid of cos i and m2 values. In each timing run, 
relativistic parameters appropriate to the given cos i and m2 were calculated and 



^ http: / /pulsar. princeton.edu / tempo 
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Table 1. Pulsar Parameters and Post-Keplerian Detections 



Pulsar 


Pb 


e 


fi 


Post-Keplerian Parameters 




(Days) 




(Mo) 


u) Shapiro 


Pb i 


J0621+1002 


8.32 


0.002 457 


0.0270 


^/ limit 




J0751+1807 


0.26 


0.000 003 


0.0010 


maybe? 


soon! 


J1713+0747 


67.83 


0.000 075 


0.0079 


V 


V 


B1855+09 


12.33 


0.000 022 


0.0056 


V 




J2019+2425 


76.51 


0.000111 


0.0107 


limit 


V 



held constant, while all other parameters were allowed to vary. The gray regions 
in figures 1 and 2 show combinations of cos i and m2 that gave acceptable fits 
at the 68% and 95% confidence levels (inner and outer contours) .□ 

For each pulsar (except J2019-I-2425), the grids of timing runs were used to 
calculate the probability distribution function (pdf) for the underlying value of 
mi. Briefly, the Ax^ value from each grid point was translated into a probability 
and, after suitable normalization, the probabilities associated with a given range 
of nil were summed, resulting in a pdf for mi. The limits on mi from each pdf 
are given in table 2. See Splaver et al. 2002 for more details on this procedure. 

3.1. PSRJ0621+1002 

The relatively high eccentricity of this system (still only e = 0.0025!) is sufficient 
to break the covariance between Pf, and uj in the timing fit. We measure uj = 
0.0116 ± 0.0008° yr"^ (la), yielding mi + m2 = 2.8 ± 0.3 Mq. This corresponds 
to the strip between dashed lines in figure la. We do not detect Shapiro delay, 
so orbits with high i are excluded, closing off the gray region towards the lower 
right portion of figure la. The low eccentricity of the system, along with the 
lack of a detectable optical companion, provides strong evidence the companion 
is a white dwarf, so m2 < 1.4 Mq. As figures la and lb show, the orbit is 
relatively close to face on, and the pulsar and companion masses are relatively 
high. Further discussion can be found in Splaver et al. 2002. 

3.2. PSRJ0751+1807 

This binary has a very short orbital period, only 6.3 hours. The relativistic 
decay rate is expected to be around Pb ~ —4 x 10~^^, depending on the system 
masses (see figure 2a for representative calculations). The current data give 
Pfe = (2 lb 8) X 10~^^, tantalizingly close to the expected value. Results of a full 
relativistic timing analysis (figure 2a) show that higher values of i are preferred, 
implying a marginal detection of Shapiro delay. The pulsar mass is not yet 
measured, though there is a strong upper limit m2 < 2.07 Mq. 

We anticipate that future observations will substantially reduce the uncer- 
tainty of Pb- (For a pulsar observed uniformly over time span t, the uncertainty 
of Pb scales as t~^-^.) This will be the first detection of relativistic decay outside 
of double-neutron-star binaries. 



^For J0621+1002, the contours are at 39% and 86% confidence levels. 
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Table 2. Mass Measurements from Timing Analysis 



Pulsar 




Piikar Mass, nii 






Best Value 


68% confidence 


95% confidence 


J0621+1002 


1.70 


1.41-2.02 


1.07-2.29 


J0751+1807 


0.55 


0.19-1.25 


0.02-2.07 


J1713+0747 


1.85 


1.49-2.26 


1.23-2.83 


B1855+09 


1.57 


1.46-1.69 


1.37-1.82 



3.3. PSRJ1713+0747 

This pulsar in a wide binary has a strong, sharp pulse, allowing us to achieve tim- 
ing precision of 200—400 ns in a single observation. Camilo, Foster, & Wolszczan 
(1994) first reported the detection of Shapiro delay. Our data give substantially 
improved constraints on the system masses and inclination (figure 2b). The 
timing data still allow a wide range of pulsar masses (table 2), but the overlap 
between parameters allowed by the timing data and the range of m2 predicted 
by the Ph—m2 relation of Tauris & Savonijc (1999) gives a tighter constraint on 
nil, 1.46— 1.61 Mq (figure 2b). (The binary evolution tracks of Podsiadlowski et 
al. 2002 give very similar predictions.) 

3.4. PSRB1855+09 

Our new data only modestly augments that of Kaspi, Taylor, &; Ryba (1994) 
and earlier work. Shapiro delay gives tight constraints on m2 and i (figure 2c), 
which are in excellent agreement with the Pb—m2 relation. 

3.5. PSRJ2019+2425 

No relativistic effects are detected. The observed x/x = 1.3 ± 0.2 x 10"-^^ s"-*^ 
constrains i < 72°. Shapiro delay is not detected, also constraining the system 
away from high i (figure 2d). Little can be said about the system masses from 
timing alone. However, if m2 is within the range allowed the Pi, — m2 relation 
(horizontal strip in figure 2d), the neutron star mass must be nii < 1.51 M©. 
Further analysis of this system is given in Nice, Splaver, & Stairs 2001. 

4. Conclusion 

Table 2 shows substantial constraints on neutron star mass nii from timing alone 
for three of the five pulsars, J0621+1002, J1713+0747, and B1855+09. In each 
case, the 68% confidence limits on mi are above the canonical neutron star mass 
of 1.35 Mq, although the 95% ranges all encompass this value. PSR J0437— 4715 
is similar, with mi = 1.58 ± 0.18 Mq (van Straten et al. 2001). The evidence is 
mounting that some pulsars are heavier than 1.35 M©, though this cannot yet 
be considered conclusive. The P5 — m2 relation is in good agreement with the 
timing data and with the relatively high pulsar masses. 

The most exciting prospect for the future is PSR J0751+1807. Although 
its mass is only weakly constrained, especially on the lower end, great advances 
will be made as the decay of this pulsar's orbit is measured in coming years. 
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Figure 1. Constraints on PSR J0621+1002 parameters. See text. 



Acknowledgments. Portions of this woric are collaborations with Zaven 
Arzoumanian, Don Baclser, Fernando Camilo, Michael Kramer, Andrea Lom- 
men, Oliver Lohmer, and Andrew Lyne. We are indebted to Joe Taylor and 
Steve Thorsett for development work on the Princeton Mark IV system and to 
Kiriaki Xilouris and Dune Lorimer for efforts in advancing pulsar astronomy 
on the post-upgrade Arecibo telescope. The Arecibo Observatory is operated 
by Cornell University under a cooperative agreement with the National Science 
Foundation. Pulsar research at Princeton University is supported by the NSF. 
IHS is supported by an NSERC UFA and Discovery Grant. 



References 



Camilo, F., Foster, R. S., & Wolszczan, A. 1994, ApJ, 437, 39 
Kaspi, V. M., Taylor, J. H., & Ryba, M. F. 1994, ApJ, 428, 713 
Kopeikin, S. M. 1996, ApJ, 467, L93 

Nice, D. J., Splaver, E. M., & Stairs, I. H. 2001, ApJ, 549, 516 
Podsiadlowski, P., Rappaport, S., & Pfahl, E. D. 2002 ApJ, 565, 1107 
Rappaport, S., Podsiadlowski, P., Joss, P. C, DiStefano, R., & Han, Z. 1995, 
MNRAS, 273, 731 

Splaver, E. M., Nice, D. J., Arzoumanian, Z., Camilo, F., Lyne, A. C, &; Stairs, 

1. H. 2002, ApJ, in press, |astro-ph/0208281 
Stairs, I. H., Splaver, E. M., Thorsett, S. E., Nice, D. J. & Taylor, J. H. 2000, 

MNRAS, 314, 459 
Tauris, T. M. & Savonije, G. J. 1999, A&A, 350, 928 
Thorsett, S. E. & Chakrabarty, D. 1999, ApJ, 512, 288 

van Straten, W., Bailes, M., Britton, M., Kulkarni, S. R., Anderson, S. B., 
Manchester, R. N., & Sarkissian, J. 2001, Nature, 412, 158 



6 



Nice, Splaver, & Stairs 



(a) J0751+1807 
0.5 




90° 80° 70° 60° 50° 40° 20°0° 
Orbit Inclination, i 



(c) B1855 + 09 
0.5 



a 

o 

'S 
to 
a 

e 

o 
u 



0.4 
0.3 
0.2 
0.1 



0.0 



_ 1 ; 1 1 1 1 ; 1 1 
m,(M,) 






; 2,00 






1.50-- 


: 

l.Off 






■ , , , 1 , , , 1 , 


. 1 , . . 1 





0.0 0.2 0.4 0.6 0.8 1.0 
cos i 

I , \ , I , I . I I I . I . 1 .11 

90° 80° 70° 60° 50° 40° 20°0'' 
Orbit Inclination, i 



(b) J1713+0747 




0.6 



to 

a 
o 

"c 
a 
ft 

E 
o 
o 



0.4 



0.2 



0.0 





1 1 1 1 1 1 1 1 1 t 1 , ,1 

/ / 1 ; ; - 

/' ^' '' /' ' 






._.3.a- 






---2.0- 




.-■1.5- 




"l.Q- 




Excluded 


< < 1 . < 



0.0 0.2 0.4 0.6 0.8 1.0 

cos i 

I , I ; I ; I I . I I I I I ill 

90° 80" 70° 60° 50° 40° 20°0° 
Orbit Inclination, i 



(d) J2019+2425 
0. 



0.6 

ta 
S 

a 0.4 



c 

as 
SX 

e 

O 
O 



0.2 



0.0 



1 1 1 1 
"< 


1/ V ' 

^/ ; 


Excluded by 


/ ///' - 


- Shapiro Delay 


















Excluded by x 

'<. . .'.> 






. 1 . , . 1 . , , 1 . < . ~ 



0.0 0.2 0.4 0.6 0.8 1.0 
cos i 

I , I , I , \ . I . I ■ I 

90° 80° 70° 60° 50° 40° 20°0° 
Orbit Inclination, i 



Figure 2. Constraints on cosi and m2. Dashed lines show repre- 
sentative values of mi. Gray regions denote the area allowed by the 
data. When present, a vertical line indicates the largest i allowed by 
a measured i, and horizontal lines indicate the range of m2 predicted 
by the Pb-m2 relation, (a) PSR J0751+1807. Solid lines show rep- 
resentative values of A- (b) PSR J1713+0747. (c) PSR B1855+09. 
(d) PSR J2019+2425. The curved solid line delimits the region in the 
upper left excluded by the non-detection of Shapiro delay. 



